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Рассмотрены вопросы специфики тепловых процессов в асинхронном двигателе при его питании несиммет- 
ричной системой напряжений. Приведены результаты компьютерного моделирования динамики нагрева. 
Определены условия возникновения опасных в тепловом отношении режимов работы. 
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Введение. Актуальной является задача моделирования процессов нагрева элементов конструк- 
ции асинхронного двигателя (АД) при его функционировании в составе вспомогательного элек- 
тропривода (ВЭП) электровоза, когда питание АД осуществляется от несимметричной системы 
3-фазных напряжений [1]. По данным статистики ОАО «РЖД» за 2010 г. выплавления обмотки 
ротора составили 33% отказов АД типа НВА-55, 25% - АНЭ-225 и 15% - АЭ 92-4. Пример конден- 
саторной схемы питания для ВЭП с АД приведен в [2]. 

Характерные особенности разработанной тепловой модели. В качестве исходных данных 
для тепловых расчетов используются результаты моделирования токов (как причин нагрева) и 
потерь в ранее разработанной в САПР ОгСАБ компьютерной модели электромеханических процес- 
сов в АД, которая позволяет рассчитывать динамические режимы АД при учете несимметрии как 
по конструкции, так и по питанию [3]. Для построения тепловой модели выбран метод тепловых 
схем с использованием принципа приближенного расчета двухмерного поля Р. Зодерберга [4]. 
Сведения о моделировании нестационарных тепловых процессов в электрических машинах, ин- 
формация по способам расчета параметров тепловых схем и справочные данные о тепловых 
свойствах электротехнических материалов приведены в [4, 5]. В тепловой модели использовано 
около 300 различных элементов типа тепловых потоков, теплоемкостей и тепловых сопротивле- 
ний. Фрагмент эквивалентной тепловой схемы пазовой части статора АД типа НВА-55 представ- 
лен на рис. 1. 

При рассмотрении динамических тепловых процессов, неодинаково развивающихся в раз- 
ных фазах, возникла необходимость учета в модели уравнительных тепловых потоков между 
фазными участками магнитопроводов, зубцовых зон, станины. Так как ротор вращается относи- 
тельно статора, то мимо различных в тепловом отношении участков обмотки статора проходит 
вся поверхность ротора (исключая работу в режиме короткого замыкания, когда ротор неподви- 
жен относительно статора). Так как ротор вращается относительно магнитного поля, то контуры 
тока, наводимые от различных фаз статора, в обмотке ротора, не привязаны жестко к частям по- 
верхности ротора. С одной стороны, если фазы статора ввиду несимметрии системы напряжений 
нагружены разными по величине токами, то в роторе из-за его вращения нет постоянной привяз- 
ки какой-либо группы стержней к определенной фазе, по этой причине не должно быть и посто- 
янно локализованного перегрева какой-либо фазы из-за несимметричного питания. С другой сто- 
роны, дефекты обмотки ротора (обрывы, раковины стержней, вырывы участков короткозамыкаю- 
щих колец и др.) привязаны именно к определенным участкам поверхности ротора, фазам. По- 
этому целесообразно пофазное рассмотрение зон в тепловой схеме ротора. 


1227 


Технические науки 













о 
Часть станины над 
частью ярма магнито- } 
провода статора, о — 
относящейся к фазе А 


К вуРН 


2 


|ВРИ 

























1 
| 
1 
| С Арн 
Часть ярма магнито- и 
провода статора, | т 
о ь 
относящаяся к фазе А а Я а ее 
1ТАРН 2 квт" 
— — сжоы 
г рн | _В о Вл Параметр, 
ии Крурн пы дамя зде и 
Пазовая часть фазы [я то лей 
«меди» статора а ой фаз 
с анизотропными у 
тепловыми 
свойствами по осям 
Часть зубцовой зоны | 
статора (7. пазов и Хх 
15и РН, 
зубцов), относящаяся 7 


к фазе А 





Е Язи? 


3 
НЕС 
оС бы 


Рис. 1. Фрагмент эквивалентной тепловой схемы пазовой части статора АД: 
Р - потери (тепловые потоки), Вт; ® - тепловое сопротивление, °С/Вт; С - теплоемкость, Дж/°С; 9; -— температура 
для различных зон (ярмо, станина, зубец, паз, клин, изоляция, металл обмотки, воздух внутри АД, воздух снаружи АДУ), °С 
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Основные потери в стали ротора в установившемся режиме обычно не учитываются ввиду 
их малости, но в процессе пуска, других переходных процессах, связанных с работой при больших 
скольжениях, такое допущение недопустимо грубо. Так как потери в стали ротора можно пред- 
ставить как выделяющиеся при протекании тока через активное сопротивление, подключенное в 
Т-образной схеме замещения параллельно традиционному контуру намагничивания [6], который, 
в свою очередь, состоит из параллельно включенных активного и индуктивного сопротивлений, 
то можно поставить основные потери в стали в зубцах ротора и в ярме ротора в соответствие ос- 
новным потерям в стали статора на основании соотношений установившегося режима короткого 
замыкания. Рассмотрим это на примере. Обозначим основные потери в стали в зубцах ротора как 
Р,„, . Их мгновенные значения рассчитываются в модели АД для каждой фазы из величины зна- 


чений потерь в стали каждой фазы, которые можно определить согласно [5] 


В 
1 п 
РН = т Ро | — п Ри - К Ве Тр> | ! (1) 


1 
где р, ‚о Удельные потери в стали 2212 толщиной 0,5 мм при индукции 1 Тл и частоте перемаг- 
ничивания 50 Гц, р, =2,2 Вт/кг; В — показатель степени, учитывающий зависимость по- 
терь в стали от частоты перемагничивания, В =1,3-—1,5 , принимаем В =1,4; п, - синхронная 
частота вращения для НВА-55, п, =1500 об/мин. 


Тогда потери в стали ротора на каждую фазу между ярмом Р,„, и зубцами Р,„„, разде- 
лятся в пропорции 
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Р, ден =: К В.Т, (2) 
Рен К т, 


где К — коэффициенты, учитывающие влияние на потери в стали неравномерности распределе- 
ния потока по сечениям участков магнитопровода и технологических факторов, К,„, =1,6, 


К, =1,8 согласно [5] для машин мощностью менее 250 квВТ. 
Для  НВА-55 индукция в ярме ротора В.› = 0,906 Тл а индукция 


=1,661 Тл, тогда Глен. _ 1,6.0,906°.51,6 














в среднем сечении зубца ротора В АО, т. е. 
м же Р„ 1,8.1,661? .16,4 
1,4 1,4 
И Рен =} г 22 72 (1,6.0,906°.51,6-+1,8.1,661°.16,4) = 72,948 | 
3 се. 1500 т 
Известно, что потери в стали НВА-55 составляют РА. = 982 Вт. При симметричном питании 
Р 
на одну фазу приходится Р„„ = т - > = 327,333 (Вт). 


В режиме короткого замыкания 








1,4 
Р,„„ =72,948 `В = 72,948 (Вт); Р„„ = 327,333 Вт; 
1500 
72,948 
Е, рн — 327 333 пен = =0, 223Р: рн . 
Обобщая последнее соотношение на другие режимы, получаем 
1,4 
П 

Р. вн = — 0, 223Р рн | о ею : (3) 

При несимметричной системе питающих напряжений 
Р —0 223 Рена ый Ё ернв с РР еенс =. п, ы , (4) 

о 3 1500 


Аналогично рассчитываются основные потери в стали ярма ротора. 

Для любого режима работы, исключая режим короткого замыкания, основные потери в 
стали в зубцах ротора и в ярме ротора разделяются в равных долях на каждую фазу. Тепловые 
потоки от поверхностных потерь в стали зубцов ротора, а также пульсационных потерь в зубцах 
ротора относятся в равных долях на каждую фазу. 

Потери в проводниках ротора в модели рассчитываются из следующих соображений. По- 
тери в проводнике пропорциональны квадрату тока и активному сопротивлению. Так как для 6е- 
личьей клетки пазовая часть одной фазы физически состоит из одного стержня, то током фазы 
будет ток через стержень. Полагаем, что при приведении обмотки ротора к трехфазной обмотке 
статора ток фазы ротора Г, пропорционален току стержня Г,... Одна лобовая часть одной фазы 


Баг 
беличьей клетки физически состоит из участка короткозамыкающего кольца между соседними 
стержнями. Каждая фаза беличьей клетки может быть описана как конструкция, содержащая 
один участок короткозамыкающего кольца между соседними стержнями, один стержень и еще 
один участок короткозамыкающего кольца между соседними стержнями. 

Ток короткозамыкающего кольца [5] 





Т ВЕ в = Ен И 
ыы 25П пр. 29п — ла 0,329 ' (5) 
7. 38 


где р - число пар полюсов НВА-55, р =2; 2, -— число зубцов ротора НВА-55, 2, = 38. 
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Активное сопротивление одного стержня 
п 
Ук5,, 16.105 .118,3.10° 
где у - удельная электрическая проводимость алюминиевого сплава АХЖ беличьей клетки ротора 


-1,26796.10* (Ом), (6) 


НВА-55 при температуре 20 °С, у = 16.10° См/м; 5,„- площадь поперечного сечения 
стержня НВА-55, 5, =118,3.10° м. 


Активное сопротивление двух участков короткозамыкающего кольца между двумя сосед- 
ними стержнями 


2%(0. и — Пт) _ 21(0,2474 - 0,045) 


И = = 0,02905.10“ (Ом). 7 
"бр мь  381610°2720.105 и (7) 

Потери в стержне можно вычислить как 
Ра =* А > 1, 26796. 10" Те (Вт), (8) 


потери в двух участках короткозамыкающего кольца между двумя соседними стержнями 
2 





1 а 1 а — га 
Ри = ми то = |: о 0,02905 .10— = 0,0883.10*12 (Вт). (9) 
558. 


Доля потерь в стержне на фазу беличьей клетки 


40-4 72 
200% - 93,5%. (10) 
10 (1,26796 +0,0883)1 


Враг 


бы 100% = 
Ра, т Р. ттдо 
Доля потерь в кольцах на фазу беличьей клетки 


0,0883 .10— 1? 


Р. 
2 1009 — Ва 100% = 6,5%. (11) 


аа 10" (1,26796 + 0,0883) 1 г 
Потери в обмотке ротора на фазу 
РоаЕе Е ееАР о 
Р. ‚= Е (12) 


2Баг 2ппдх2 
Приведенные выше соотношения справедливы при симметрии фаз обмотки ротора (без 
дефектов стержней). В случае, когда в одной из фаз, например А, оборваны п, стержней, поте- 


ри в обмотке ротора на фазу могут быть рассчитаны для фазы А 





р _ Ра +Рив + Риме Ри | _ ты 13) 
25РНА ^^ г 
т О т о 2, 
Аналогично для фазы В: 
р Род Рив Рае Ры | | пт ва 
25РНВ ь 
т Ра т Р, пох 74 
Если в фазе С нет оборванных стержней, то потери в ее пазовой части будут 
р 
а 2Ь 
Р.рнс = (Ре + Р.рнв + Ррнс ) р р — Р.брив — Роба. (15) 
2Баг т 2пп9х2 


Если обрывы стержней есть лишь в одной фазе, например А, то для фаз без разрывов 
стержней: 
р 


1 а 
Р.рнс = 5 (Ро + Р.рнв + РЬрнс — Робера. (16) 
2раг 


2ипдх2 
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Если в каждой фазе есть разрывы стержней, то используются выражения (например, ко- 
гда оборвано минимальное количество стержней п,. в фазе ©): 


Р 
= 2Ба ь 
Р, рис я” Вы ТЕ Р. прив те Ра ) р Г. = РО РЕВ — Ре Г (17) 
2Баг к 2ппдх2 

ро феи + Рав Ре _Ры  ||_ _ПыТ | а 
25РНА — р р р Рраеищет | (18) 

2раг 2 пп х2 р) оС 

р Рем +В Ре Гы |1 П,вТ 
25РНВ — т р р 7 РИ я (19) 

2аг г 2ипдх2 с 


Чтобы чрезмерно не усложнять тепловую модель и не замедлять ход решения на ЭВМ, 

приняты некоторые допущения. Например, тепловое сопротивление между внутренним воздухом 
машины и наружным воздухом, определяемое подогревом внутреннего воздуха, принято постоян- 
ным, хотя фактически оно является переменным и зависит от частоты вращения ротора (скорости 
перемешивания воздуха). 
Определение адекватности тепловой модели. Вычисления направлены на установление 
адекватности тепловой модели экспериментальным результатам. С учетом реальной длительности 
(3,540 ч) установления тепловых процессов в НВА-55 моделирование таких процессов в 
ОГСАР 9.2 от нулевого перегрева до установившегося представляется затруднительным: требует- 
ся машинное время порядка нескольких суток, значительный ресурс оперативной памяти ЭВМ, 
предельный размер файла результатов моделирования - 2 Гбайт. Поэтому проведены вычисли- 
тельные эксперименты, позволяющие оценить адекватность модели при существенно меньших 
затратах машинного времени: 

— моделирование адиабатического процесса нагрева обмоток (стоянка под током с горяче- 
го состояния) при симметричном синусоидальном питании напряжением 380 В с оценкой средней 
температуры обмотки статора по истечении 18 с; 

— моделирование длительной работы АД для привода вентилятора, питаемого по конден- 
саторной схеме, с оценкой промежуточных результатов нагрева участков станины, расположен- 
ных над железом статора. 

В качестве граничных условий в начале моделирования процесса адиабатического нагре- 
ва для имитации горячего состояния машины приняты данные, базирующиеся на результатах ис- 
пытаний НВА-55 ОАО «ВЭлНИИ» при номинальной нагрузке (табл. 1). 

Таблица 1 

Граничные условия для моделирования в начале процесса адиабатического нагрева, 

имитирующие горячее состояние НВА-55, базирующиеся на результатах испытаний НВА-55 
при номинальной нагрузке 


И: и, в Температура обмотки Температура обмотки Температура Температура под- 
- - - статора, °С ротора, °С корпуса, °С шипников, °С 





381 


Принятые граничные условия приведены в табл. 2. Сравнение результатов моделирования 
процесса адиабатического нагрева НВА-55 (стоянка под током с горячего состояния) с опытными 
данными при температуре окружающего воздуха 20 °С дано в табл. 3. Наблюдается удовлетвори- 
тельное соответствие результатов моделирования данным опыта. Согласно [7], перегрев алюми- 
ния беличьей клетки ротора не должен быть выше 300 °С, т. е. температура не должна превы- 
шать 350 °С, что в сопоставлении с данными табл. 3 косвенно также свидетельствует об адекват- 
ности модели. Данное температурное ограничение объясняется значительным увеличением пла- 
стичности алюминиевых сплавов при перегреве выше 300 °С. 
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Таблица 2 
Принятые граничные условия моделирования в начале процесса адиабатического нагрева 
для имитации горячего состояния НВА-55 


Температура зубцов Температура ярма Температура зубцов Температура ярма ротора, вала, подшипниковых 


статора, °С статора, °С ротора, °С щитов, внутреннего воздуха, °С 


50 


Экспериментальной базой для определения адекватности моделирования длительной ра- 
боты АД для привода вентилятора (мотор-вентилятора), питаемого по конденсаторной схеме, с 
оценкой промежуточных результатов нагрева станины над железом статора стали данные испы- 
таний ОАО «ВЭлНИИ» на электровозе ЗЭС5К-047, а именно кривые изменения перегрева самой 
нагретой части станины НВА-55 над температурой окружающего воздуха во времени. Ввиду 
большой длительности процесса нагрева (до 4 ч) и ограниченности вычислительных ресурсов 
(моделирование проводилось на ЭВМ с процессором типа РепНит 4 и оперативной памятью 2 
Гбайт при предельном размере файла результатов моделирования 2 Гбайт) для сравнения дос- 
тупны начальные участки кривых нагрева (см. рис. 2). 





Таблица 3 
Сравнение результатов моделирования процесса адиабатического нагрева НВА-55 
(стоянка под током с горячего состояния) 
с опытными данными при температуре окружающего воздуха 20 °С 


Температура по- Температура корот- Температура Температура лобо- Температура пазо- 
верхности корпуса козамыкающих колец стержней обмотки вых частей обмотки | вых частей обмотки 
над железом стато- обмотки ротора, °С ротора, °С статора, °С статора, °С 

ра, °С (модель) 
(модель) (модель) (модель) (модель) 


Длительность 
процесса, с 


363,4 343,6 166,6 158,8 


Средняя температура обмотки ротора, °С Средняя температура обмотки статора, °С 
(модель) 


модель ОПЫТ 





353,5 


“ 


и нагретои части корпуса 


“ 
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Рис. 2. Сравнение опытных и расчетных данных о нагреве станины НВА-55 
при питании по конденсаторной схеме несимметричной системой напряжений 
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Из рис. 2 видно, что процесс нагрева на модели идет несколько быстрее, чем по опытным 
данным (средняя относительная ошибка 20%). Различие может быть скорректировано за счет 
увеличения теплоемкости внутреннего воздуха машины путем увеличения его объема (при моде- 
лировании не учтен обмен воздухом с окружающей средой). В пользу именно такого объяснения 
природы расхождения результата опыта и моделирования на рис. 2 говорит то обстоятельство, 
что при моделировании адиабатического нагрева, где за время процесса практически не успевает 
происходить обмен внутреннего воздуха машины с наружным, расхождение расчетной и опытной 
величин среднего превышения температуры обмотки статора над температурой окружающего 
воздуха меньше и имеет противоположный знак (относительная ошибка составляет минус 11%). 
Отметим, что тепловая модель станины НВА-55 (см. рис. 1) носит приближенный характер ввиду 
довольно сложных путей теплопередачи от магнитопровода статора к корпусу, что обусловлено 
спецификой конструкции этого АД. 

Результат моделирования тепловых процессов при включении одного мотор-вентилятора 
в длительном режиме по конденсаторной схеме питания показан на рис. 3. 


Тпер, °С 


п, об/мин 


50 





Рис. 3. Результат моделирования тепловых процессов при включении одного АД привода вентилятора 
в длительном режиме по конденсаторной схеме (графики перегревов самой нагретой фазы А): 
1- короткозамыкающие кольца обмотки ротора; 2- стержни обмотки ротора; 3 -— лобовые части обмотки статора; 
4- пазовая часть обмотки статора; 5-— поверхность станины; 6 - частота вращения ротора //100, об/мин 


Известно [8, 9], что частота включений АД привода компрессора (мотор-компрессора) мо- 
жет достигать 40—60 в час, что сопровождается значительными токами и, как следствие, выделе- 
нием тепла. Результаты моделирования режима ПВ мотор-компрессора, питаемого по конденса- 
торной схеме, при емкости конденсаторов 70% от номинальной (это минимальный предел сниже- 
ния емкости от номинала, допускающий, согласно результатам моделирования, пуск мотор- 
компрессора при напряжении на обмотке собственных нужд тягового трансформатора не ниже 
380 В) показаны на рис. 4 (ПВ 50% при длительности цикла 8 с). В режиме работы, показанном на 
рис. 4, перегрев алюминия беличьей клетки ротора до 300 °С может быть достигнут за время по- 
рядка 500—600 с. При тяжелых затяжных пусках ввиду большей продолжительности протекания 
пусковых токов это время сократится. 


1233 


Технические науки 


Тпер, °С 


п, об/мин 


200 


100 





Рис. 4. Результаты моделирования режима ПВ мотор-компрессора, питаемого по конденсаторной схеме 
при емкости конденсаторов 70% от номинальной (графики перегревов самой нагретой фазы А) 


Определение опасных в тепловом отношении режимов работы АД. Поставлена задача 
определить, при каких условиях в случае питания АД типа НВА-55 по конденсаторной схеме при 
номинальной емкости конденсаторов за относительно короткий промежуток времени достигается 
перегрев алюминиевого сплава марки АХЖ беличьей клетки ротора до 300 °С. С помощью разра- 
ботанной тепловой модели эти условия выявлены. Индикатором появления состояния обмотки 
ротора, которое влечет за собой весьма быстрое (за 8-20 с) нарастание перегрева участков бе- 
личьей клетки из сплава АХЖ до опасных значений 300 °С и более, является возникновение при 
питании по конденсаторной схеме небаланса междуфазных напряжений на частоте основной гар- 
моники с величиной небаланса 360-380 В при установившемся режиме работы привода. Опасные 
по нагреву режимы сопровождаются весьма значительными величинами токов фаз АД, причем в 
одной из фаз ток в 7-9 раз превосходит номинальный. Коэффициент небаланса междуфазных 
напряжений определялся согласно [10] 


Кебмь = 5 100%, (20) 


Н 


где И’, Ц. - наибольшее и наименьшее действующие значения из трех междуфазных напряже- 
ний; И — номинальное значение линейного напряжения АД, Ш. = 380В. 


Как показало компьютерное моделирование пуска до установившегося режима единично- 
го АД, такие опасные режимы работы складываются в результате обрыва 21-33% стержней бе- 
личьей клетки при определенных вариантах их распределения по фазам ротора. Они характери- 
зуются возможностью достижения температуры плавления алюминиевого сплава в отдельных 
участках ротора за время 20—40 с (температура плавления технического алюминия 658 °С) при 
перегреве поверхности корпуса АД 0,1-0,3 °С. При этом АД работает в установившемся режиме 
со скольжением 7—36% (скольжение НВА-55 в номинальном режиме составляет 3,9%). Коэффи- 
циент небаланса междуфазных напряжений составляет Кв = 96-100%. Для развития опасного 


режима достаточно повреждения хотя бы одного стержня дополнительно к ранее оборвавшимся, 
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причем до этого критического обрыва коэффициент небаланса междуфазных напряжений был в 
пределах К..„, = 24-36%. Для случая питания единичного НВА-55 по конденсаторной схеме при 


отсутствии дефектов обмотки ротора К„„„„ = 23,7%. 


Пример множественного повреждения стержней беличьей клетки НВА-55 показан на рис.5 
(фотография ротора АД, эксплуатировавшегося на электровозе типа ЭП1М). Можно рассмотреть 
выплавления и разрывы стержней (не менее 12 штук) почти во всех пазах, которые видны на фото. 





Рис. 5. Ротор АД НВА-55 с множественными повреждениями стержней беличьей клетки, 
эксплуатировавшийся на электровозе типа ЭП1М с конденсаторной системой питания АД 


Среди множества причин возникновения и развития дефектов обмотки ротора из алюми- 

ния и его сплавов называют остаточные механические напряжения в процессе изготовления бе- 
личьей клетки, что подробно описано в [11]: так как коэффициент температурного расширения 
алюминия в 2 раза выше, чем стали, при остывании алюминиевой заливки сердечника возникают 
механические силы в стержнях беличьей клетки, которые могут превосходить по величине предел 
текучести алюминия, приводя к податливости стержней. 
Заключение. Разработанная тепловая модель АД является адекватной и может быть использо- 
вана для исследования тепловых процессов, включая случаи несимметричного питания и конст- 
руктивной несимметрии АД, в том числе для выявления опасных в тепловом отношении режимов 
работы электропривода. Например, режим ПВ при частых включениях может быть значительно 
более напряженным в тепловом отношении, чем длительный режим. Индикатором опасных по 
нагреву режимов, вызванных дефектами стержней беличьей клетки ротора, является коэффици- 
ент небаланса междуфазных напряжений АД, который целесообразно периодически контролиро- 
вать в процессе эксплуатации по показаниям бортовой микропроцессорной системы локомотива. 
Критическое значение этого коэффициента может быть установлено по результатам сопоставле- 
ния измерений в процессе эксплуатации со статистикой отказов вспомогательных АД. 
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